Journal  of  Power  Sources  196  (2011)  692-699 


ELSEVIER 


Contents  lists  available  at  ScienceDirect 

Journal  of  Power  Sources 

journal  homepage:  www.elsevier.com/locate/jpowsour 


in 

SllbudtS 


Tortuosity  in  anode-supported  proton  conductive  solid  oxide  fuel  cell  found  from 
current  flow  rates  and  dusty-gas  model 

Chih-Long  Tsai*,  V.  Hugo  Schmidt 

Department  of  Physics,  Montana  State  University,  Bozeman,  MT  5971 7,  USA 


ARTICLE  INFO 


ABSTRACT 


Article  history: 

Received  22  June  2010 
Received  in  revised  form  30  July  2010 
Accepted  3  August  2010 
Available  online  8  August  2010 


Keywords: 
Proton-conductive 
Solid  oxide  fuel  cell 
Tortuosity 
Dusty-gas  model 


Ba(Ceo.8Yo.2)03_5  anode-supported  proton  conductive  solid  oxide  fuel  cells  were  fabricated  and  tested. 
By  changing  the  H2  partial  pressure  at  the  anode  side,  the  effect  of  anodic  concentration  polarization  on 
open-circuit  voltage  of  the  cell  was  observed.  Saturation  current  densities  under  concentration  polar¬ 
ization  were  obtained  from  different  anode  thickness  cells  and  were  used  for  tortuosity  calculation.  The 
calculation  is  based  on  the  dusty-gas  model  which  includes  Knudsen  diffusion  and  Stefan-Maxwell  equa¬ 
tion  terms.  The  tortuosity  value  for  our  supporting  anode  is  1.55  ±  0.1  which  is  in  a  physically  reasonable 
range  for  modern  porous  anode  materials.  The  tortuosity  that  we  found  is  independent  of  the  cell  testing 
temperature  and  anode  thickness,  which  is  consistent  with  the  fact  that  tortuosity  is  a  geometric  factor  of 
the  anode  structure.  The  derived  equation  also  can  be  used  for  predicting  the  effect  of  varying  the  anode 
thickness,  porosity  and  pore  size.  Also,  the  concentration  of  the  gases  as  a  function  of  position  across  the 
anode  is  determined. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Strong  efforts  of  modern  solid  oxide  fuel  cell  (SOFCs)  designs 
are  in  progress  to  extend  the  operating  temperature  down  to  the 
600-800  °C  range.  The  low  operating  temperature  would  increase 
the  benefit  of  using  low  cost  interconnect  material,  such  as  stainless 
steel,  and  also  reduce  problems  associated  with  thermal  expan¬ 
sion,  atomic  migration  and  corrosion.  Accordingly,  the  electrolyte 
is  made  as  thin  as  possible  to  reduce  ohmic  loss  from  transferring 
ions  through  electrolyte.  The  electrolyte  thickness  is  usually  down 
to  20  pim  range  and  requires  that  an  electrode  layer,  usually  the 
anode,  must  be  made  thick  enough  to  mechanically  support  the 
cell. 

With  a  thicker  anode,  gas  concentrations  at  the  anode- 
electrolyte  interface  need  to  be  calculated  carefully  because  they 
strongly  affect  the  terminal  voltage  V  as  a  function  of  electrolyte 
current  density  i,  especially  in  the  high  current  density  range.  A  key 
parameter  in  determining  the  pressure  gradient  is  the  tortuosity  r. 
The  definition  of  ‘tortuosity’  is 

typical  diffusion  path  between  two  points  through  the  pores 

tortuosity  =  — - — — f- - - - -A- - - - - 

straight  distance  betwee  the  same  two  points 
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and  sometime  is  confused  with  ‘tortuosity  factor’  which  is  defined 
as 

tortuosity  factor  =  (tortuosity)2 

Measurement  of  its  value  is  needed  for  determining  the  quality  of 
the  anode  pore  configuration  and  necessary  for  any  V(i)  modeling. 
Many  modern  SOFC  models  do  not  adequately  calculate  the  value  of 
r  and  invoke  anode  tortuosities  in  the  range  as  high  as  10-17  [1,2]. 
Such  high  tortuosities  do  not  seem  physically  reasonable  when  the 
porosity  is  usually  higher  than  30%  of  the  supporting  anode. 

The  impact  of  such  high  tortuosity  (10-17)  is  to  produce  con¬ 
centration  polarization  easily  when  running  the  cell.  With  this 
misleading  analysis,  one  might  think  sacrificing  the  structural 
integrity  to  decreasing  the  anode  thickness  or  increasing  its  poros¬ 
ity  is  advantageous.  Flowever,  recent  studies  of  SOFC  supporting 
anodes  using  focused  ion  beam  scanning  electron  microscopy 
(FIB-SEM),  X-ray  computed  tomography  and  gas  counter-diffusion 
provide  evidence  that  the  tortuosity  for  typical  supporting  anodes 
is  1.33-4.0  [3-5].  The  present  work  provides  another  way  for 
finding  tortuosity  of  the  supporting  anode  for  proton  conductive 
solid  oxide  fuel  cells  (FI-SOFCs)  by  finding  saturation  current  den¬ 
sity  under  concentration  polarization  of  H-SOFC  and  using  the 
dusty-gas  model.  Tortuosity  calculations  from  our  tested  cells  indi¬ 
cate  that  high  tortuosity  values  may  not  be  correct.  The  results 
also  provide  more  accurate  information  for  SOFC  design  and  V(i) 
modeling. 
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Nomenclature 

a 

anode  forward  attempt  current  density  (A cm-2) 

a0 

open-circuit  value  of  a  (A  cm-2) 

001 

open-circuit  value  of  a  whenpH2  =  1  atm  (A  cm-2) 

<4/3 

unit  volume  of  oxygen  ion  in  unit  lattice  of  elec¬ 
trolyte  material  (cm3) 

Dia 

Knudsen  diffusion  coefficient  for  component  i 
(cm2  s-1) 

Dij 

binary  diffusion  coefficient  (cm2  s-1 ) 

% 

binary  diffusion  coefficient  when  total  pressure 

Pi  =  1  atm  (cm2  s-1) 

fa 

fraction  of  oxygen  sites  at  the  anode-electrolyte 
interface 

i 

SOFC  current  density  (A cm-2) 

las 

saturation  current  density  under  anodic  concentra¬ 
tion  polarization  (A cm-2) 

f 

molecular  flow  densities  of  gas  i  (molec  s-1 ) 

Ji,m 

metered  inflow  rates  of  gas  i  (molec  s-1 ) 

Ji,net 

net  inflows  into  the  plenum  of  gas  i  (molec  s-1 ) 

k 

Boltzmann’s  constant,  1.38  x  10-23  (JK-1) 

Mi 

molar  mass  of  the  diffusing  gas  i  (g  mol-1 ) 

mi 

molecule  mass  of  component  i  (g  molec-1 ) 

ni 

concentration  of  component  i  (molec  cm-3) 

n\ 

total  gas  concentration  in  the  plenum  when 

P  =  1  atm  (molec cm-3) 

nlP 

concentration  of  component  i  in  the  plenum 
(molec  cm-3) 

0  H2,as 

H2  concentration  at  the  anode-electrolyte  interface 
under  saturation  current  density  (molec cm-3) 

Tli,p,as 

concentration  of  gas  i  in  the  plenum  at  the  outer 
anode  surface  under  saturation  current  density 
(molec  cm-3) 

nH,ele 

proton  concentration  in  the  electrolyte  material 

(cm-3) 

P 

total  pressure 

Q 

charge  transfer  per  reaction,  3.2  x  10-19  (C) 

f 

mean  pore  radius  (m) 

T 

absolute  temperature  (K) 

Ua 

enthalpy  released  when  the  reaction  at 
anode-electrolyte  interface  occurs  (Jkgel-1) 

V 

SOFC  terminal  voltage  (V) 

Va 

potential  difference  across  the  anode-electrolyte 
interface  (V) 

vele 

probability  of  oxygen  ion  vacancy  in  electrolyte 
material 

W 

anode  thickness  (cm) 

Greek  letters 

average  collision  diameter  of  molecules  i  and  j  (A) 

X 

tortuosity 

<t> 

porosity  of  supporting  anode 

i 2 

dimensionless  collision  integral 

2.  Experimental  procedures 


with  Cu  Ka  (X  =  1.5418  A)  was  used  for  checking  the  formation  of 
Ba(Ce0.8Yo.2)03-a-  auto-grinding  and  calcining  processes  were 
repeated  until  a  single  phase  of  the  material  was  confirmed  by  the 
X-ray  diffraction. 

To  prepare  the  supporting  anodes,  the  Ba(Ceo.8Yo.2)03_5  powder 
was  weighed  and  mixed  with  NiO  at  volume  ratio  1  to  2.  In  addi¬ 
tion  to  ceramic  powders,  8  wt%  cornstarch  of  the  total  solid  load 
was  added  to  the  powder  and  auto-grinded  for  2  h  to  serve  as  pore 
former.  The  powders  were  uniaxial  die-pressed  with  a  3/4"  die  to 
form  green  pellets  to  serve  as  the  supporting  anode.  Then,  the  anode 
interlayer,  with  a  1:1  volume  ratio  of  Ba(Ce0.8Yo.2)03_s  to  NiO,  and 
the  electrolyte  were  paint  brushed  on  the  supporting  anode  using 
the  inks  prepared  by  mixing  the  solid  powders  with  alpha  terpinol, 
ethylcellulose,  oleic  acid  and  xylene  using  a  three  roll  mill.  After  the 
drying  of  electrolyte  ink,  the  pellets  were  re-pressed  by  the  uniaxial 
die  press  and  sintered  in  air  at  a  temperature  of  1400  °C  for  5  h  to 
form  a  well-bonded  electrolyte-anode  structure.  Cathode  material, 
La0.8Sro.2Mn03  (LSM),  was  then  applied  on  the  sintered  pellets  by 
paint  brush  and  fired  to  1000  °C  in  air  for  2  h  to  form  a  good  bond 
between  cathode  and  electrolyte. 

A  home-made  seal-less  testing  system  made  from  Inconel  600 
was  used  for  the  measurements,  Fig.  1 .  Silver  mesh  and  nickel  foam 
were  used  as  cathode  and  anode  current  collectors,  respectively. 
The  total  input  gas  flow  rates  on  each  side  were  controlled  by 
MKS  mass  flow  controllers  at  200  ml  min-1  in  all  experiments.  The 
water  vapor  concentration  in  fuel  gas  was  about  3%  by  flowing  the 
fuel  gas  through  a  water  bubbler  at  room  temperature.  Measure¬ 
ments  were  carried  out  at  700  and  800  °C  in  ambient  pressure.  All 
cells  were  reduced  in  situ  at  high  temperature  in  a  60%  H2  +40% 
N2  mixture  for  more  than  an  hour  prior  to  the  measurements.  All 
performance  of  the  cells  was  measured  using  various  fuel  gas  mix¬ 
tures,  i.e.  at  various  partial  pressures  of  H2/N2  ratio  in  fuel  gas  while 
the  total  flow  rate  of  the  gas  was  kept  constant  at  200  ml  min-1 . 
The  current  densities  were  calculated  based  on  the  cathode 
area. 

Porosity  of  the  supporting  anode  was  measured  using 
Archimedes’  method.  Each  of  the  tested  cells,  after  reducing,  was 
broken  into  two  pieces.  One  of  the  pieces  was  polished  to  erase 
the  cathode  and  electrolyte  layers  and  cleaned  with  an  ultrasonic 
bath  in  ethanol.  The  sample  was  then  kept  in  a  dry  oven  at  95  °C 
for  more  than  2  h  to  evaporate  the  ethanol.  Dry  weight,  Wdry ,  wet 
weight,  Wwet,  and  weight  saturated  with  ethanol,  Wsat,  were  mea¬ 
sured  using  a  high  accuracy  balance.  Prior  to  measuring  Wwet  in 
ethanol,  the  sample  was  immersed  in  ethanol  and  kept  in  vacuum 
for  10  min  to  remove  possible  air  in  the  pores.  The  Wsat  was  mea¬ 
sured  in  air  soon  after  the  surface  of  the  sample  was  shaken  dry. 
‘Porosity’  was  calculated  using  the  equation 

Porosity  TO  *'00%  <’> 

The  other  part  of  the  cell  was  hardened  in  an  epoxy  and  pol¬ 
ished.  Field  Emission  Scanning  Electron  Microscopy  (SUPRA™  55 
Versatile  High  Performance  FE-SEM,  Zeiss)  was  used  to  examine 
the  microstructure  of  the  cell.  The  mean  pore  radius  of  the  cell  was 
determined  by  quantitative  measurements  of  the  pore  size  on  SEM 
images. 


The  Ba(Ce0.8Y0.2)O3_5  powder  was  prepared  from  solid  state  3.  Dusty-gas  model  for  gas  flow  in  anode-supported  H-SOFC 

reaction.  Stoichiometric  amounts  of  barium  carbonate  (BaC03,  pores 

99.8%),  cerium  oxide  (Ce02,  99.97%),  and  yttrium  oxide  (Y203, 

99.99%)  were  mixed  by  an  agate  auto-grinding  machine  for  two  The  multi-gas  diffusion  process  in  pore  structure  is  gen- 

hours.  Then,  the  powder  was  calcined  at  1100  °C  for  15  h  to  form  erally  described  by  the  dusty-gas  model,  which  includes  the 
the  perovskite  phase.  An  X-ray  diffraction  (Scintag,  XGEN-400)  Stefan-Maxwell  equation  and  Knudsen  terms.  The  equation  in 
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molecular  units  is  [6] 

Jj  C litij-JjtiiW  9fij 

+  DU,P,  9x  1  J 

where  J;  and  Jj  are  molecular  flow  densities  of  components  i  and 
j  (molecs-1),  respectively,  DKi  is  the  Knudsen  diffusion  coefficient 
for  component  i,  is  the  binary  diffusion  coefficient  at  Pi  =  1  atm, 
and  the  nkT- type  terms  in  the  numerator  divided  by  P\  convert  Dy! 
to  Dy ,  the  binary  diffusion  coefficient,  at  the  actual  total  pressure  at 
any  position  x  along  the  anode  pore,  nz  and  rij  are  the  concentration 
of  components  i  and  j,  k  is  Boltzmann’s  constant,  1 .38  x  1 0-23  J  K-1 , 
and  T  is  the  absolute  temperature.  It  is  important  to  clarify  that  x 
in  Eq.  (2)  is  the  coordinate  along  a  typical  gas  diffusion  path,  not 
across  the  anode.  Therefore,  the  total  diffusion  path  length  is  tw, 
the  product  of  the  tortuosity  and  the  anode  thickness. 

The  Knudsen  diffusion  coefficient  considers  the  collisions 
between  gas  molecules  and  the  wall.  The  equation  for  the  Knudsen 
diffusion  coefficient  from  the  kinetic  theory  of  gases  is 


(3) 


where  r  is  the  mean  pore  radius  and  M,-  is  the  molar  mass  of  the 
diffusing  gas.  The  calculated  Knudsen  diffusion  coefficients  for  H2, 
N2  and  H20  at  700  and  800  °C  are  listed  in  Table  1(a). 


Table  1 

Calculated  (a)  Knudsen  and  (b)  binary  diffusion  coefficients  for  various  gases  at  700 
and  800  °C  with  r  =  1 .41  (Jim. 


(a) 

DKi  (cm2/s) 

h2 

n2 

h2o 

800  °C 

31.68 

8.47 

10.56 

700  °C 

30.17 

8.06 

10.06 

(b) 

Dy!  (cm2/s) 

h2-n2 

h2-h2o 

n2-h2o 

700  °C 

5.35 

6.65 

1.87 

800  °C 

6.10 

7.70 

2.21 

The  binary  diffusion  coefficient,  Dy,  is  calculated  using  the 
Chapman-Enskog  equation  from  Cussler  [7] 


1.86  x 


Dij  = 


lO-3  T2/3  f  1,1 

1U  1  '  Mj 


\  V2 


(4) 


where  Q  is  a  dimensionless  collision  integral,  based  on  the 
Lennard-Jones  potential,  cry  is  the  average  collision  diameter  of 
molecules  i  and  j  (in  A),  M*  and  Mj  are  molar  masses  of  diffusion 
gases  i  and  j,  respectively,  and  P  the  total  pressure  (in  atm).  Using 
Table  2.2  and  2.3  in  Cussler  for  Q  and  cry  and  total  pressure  P  =  1  atm, 
the  calculated  Dyi  for  various  gases  at  temperature  700  and  800  °C 
are  listed  in  Table  1(b). 

The  experiments  were  performed  with  flow  of  H2/N2  gas  mix¬ 
tures  to  the  anode  of  Ba(Ceo.8Yo.2  )03_5  anode-supported  SOFCs.  The 
H2/N2  gas  mixture  picked  up  H20  (steam)  through  a  water  bubbler 
so  that  it  is  a  ternary  system,  or  a  binary  system  when  1 00%  H2  was 
used.  The  flow  density  of  H2  is  Jh2  =  ir2/0g,  where  i  is  the  current 
density  (A  cm-2  =  Ccnrr2  s-1)  in  the  solid  electrolyte  and  q  is  the 
charge  (3.2  x  1 0-19  C  molec-1 )  carried  per  H2  molecule  annihilated 
in  the  reaction  at  the  anode/electrolyte  interface. 

The  t2/0  factor  is  an  enhancement  factor  by  which  the  flow 
density  JHl  is  enhanced  compared  to  its  value  if  the  anode  were 
completely  porous  (i.e.  0  =  r  =  l).  To  derive  this  result,  consider 
a  simple  model  for  which  the  anode  has  pores  of  circular  cross- 
section  which  are  tilted  at  an  angle  0  away  from  the  normal  to  the 
anode-electrolyte  interface.  The  tortuosity  r  (actual  path  length 
divided  by  anode  thickness)  is  1  /cos 6.  The  flow  velocity  component 
perpendicular  to  the  interface  is  a  fraction  cos 0  of  the  velocity  in 
the  pore,  so  this  contributes  a  factor  r  to  the  enhancement  factor.  A 
plane  parallel  to  that  interface  has  pore  area  fraction  0  compared  to 
the  total  area,  thereby  contributing  a  factor  1  /0  to  the  enhancement 
factor.  The  pore  area  intercepted  by  this  plane  has  elliptical  shape, 
but  the  flow  is  perpendicular  to  the  pore’s  circular  cross-section 
which  is  smaller  by  a  factor  cos 0,  so  this  contributes  another  factor 
r  to  the  enhancement  factor.  In  earlier  work  [6]  we  overlooked  this 
final  factor  and  incorrectly  used  t/0  as  the  enhancement  factor. 
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We  know  that  there  is  no  N2  and  H2O  flow  in  the  supporting 
anode,  i.e  Jn2  =  0  and  Jh2o  =  0.  Eq.  (2)  can  then  be  written  as 


9nH 2  ( it2 \  (  1  n^kT  nHlkT  \ 

dx  \4>q  J  \A<,h2  +  Dh2,n2,iPi  +  ^h2,h2o,i^i  ) 

9nN  it2  nN  kT 

— - — —  and 
dx  (pq  DN2,h2,i^i 

9nH2o  _  h2  n^okT 
dx  (pq  E>h2,h2 0,1^1 


(5a) 

(5b) 

(5c) 


In  reviewing  mass  transport  modeling  literature,  none  of  the 
paper  apply  this  enhancement  factor  in  their  calculation  [8-11]. 
Nevertheless,  an  effective  diffusion  coefficient  was  used  for  the  cal¬ 
culations  without  further  explanation,  which  is  in  general  defined 
as 

of  =  t Dj  (6) 

where  and  D*  are  effective  and  general  diffusion  coefficient 
of  species  z,  respectively.  However,  after  carefully  reviewing  the 
FIB-SEM  and  XCT  articles  [3,5],  and  considering  discussion  of 
diffusion-based  tortuosity  analysis  by  Moldrup  et  al.  [12]  and  of 
confusion  between  tortuosity  and  tortuosity  factor  by  Epstein  [13], 
it  appears  that  the  “tortuosity”  in  these  articles  and  in  Eq.  (6)  should 
be  “tortuosity  factor”  which  is  r2.  Then  our  enhancement  factor 
r 2\(p  agrees  with  the  ones  used  in  Refs.  [3,5]  and  [8-1 1  ]. 

Adding  Eqs.  (5a)-(5c)  yields  the  equation  for  total  concentration 
n  and  its  solution 


3  n 
dx 


it * 


1 


(pq  dk,h2 


and  n  —  Tit0tai,p 


(pqD  k,h2 


—  Ktotal,p  CU2  ? 

(8) 


where  ntota/p  is  the  total  gas  concentration  in  the  plenum  and  x =yt 
is  used  in  the  equation,  where  y  is  position  on  a  straight  line  going 
across  the  anode.  From  the  ideal  gas  law,  P=nkT,  we  see  that  the 
total  pressure  P  decreases  linearly  with  x  from  the  anode-plenum 
interface  to  the  solid  electrolyte. 

We  can  solve  Eq.  (5a),  (5b)  and  (5c)  for  nH2 ,  «n2  and  rzH2o 

nN2  =nN2,pexp  (n|^N2  |)  =  nN2,p exp(cN2r3)  (8a) 


Hh20 


nH2,p  exp(cH2oT3)  and 


(8b) 


Hh2  =  riH2,p -Ch2t3  -nN2,p(exp(cN2r3)-  1) 

-nH2o,p(exp(cH2or3)- 1)  (8c) 

where  rzN2?p,  nH2o,p  and  rzH2,P  are  concentrations  for  N2,  H20  and 
H2  in  the  plenum,  respectively,  n.\  is  the  total  gas  concentration 
in  the  plenum  when  P=  1  atm  =  1.015  x  105  N  m-2,  n\  =P\lkT  and 
ni  =7.559  x  1018  cm-3  for  T=973  K  and  m  =6.855  x  1018  cm-3  for 
T  =1073  K.  We  see  that  nN2  and  nH2o  increase  exponentially  along 
the  pore  from  plenum  to  electrolyte. 

If  the  cell  current  density  is  saturated  by  concentration  polar¬ 
ization  at  the  anode,  designated  as  ias,  then  z  become  ias  in  the 
“constants”  Ch2,  Cn2  and  cH2o,  which  are  really  linear  functions  of 
position  x/r  across  the  anode.  For  x/r  =  w  (the  anode  thickness) 
at  the  anode-electrolyte  interface,  they  really  become  constants 
and  are  designated  as  CH2,as»  cn2,as  and  cH2o,as  for  the  saturated 
condition. 


To  find  Hh2  at  the  anode-electrolyte  interface,  we  use  results 
from  our  analysis  of  H-SOFC  voltage  V  as  a  function  of  cur¬ 
rent  density  z.  The  analysis  of  the  I-V  curve  for  H-SOFC  is  based 
on  the  concept  of  exchange  current  density,  z0.  The  exchange 
current  density  is  defined  for  open-circuit  operation  in  which 
the  equal  and  opposite  current  densities  in  the  system  result 
from  the  equal  forward  and  reverse  reaction  rates  at  the  inter¬ 
face  of  anode-electrolyte,  H2*>2H+  +  2e_,  cathode-electrolyte, 
2H+ +  02_  **H20,  and  cathode-pore,  l/202  +  2e_  ^>02_.  In  our 
model,  we  propose  a  two-step  cathode  reaction  mechanism 
instead  of  a  one-step  mechanism.  For  a  one-step  cathode  reac¬ 
tion,  4H+  +  02  +4e~  *>2H20,  an  02  molecule  is  required  to  react 
with  four  protons  and  four  electrons  at  a  cathode  TPB  site  to 
form  2H20  molecules.  However,  the  concentration  of  protons  in 
the  electrolyte  is  low  (<20%),  which  makes  the  probability  of  four 
protons  simultaneously  being  within  reaction  distance  is  too  low 
to  provide  the  dominant  reaction  mechanism.  Accordingly,  there 
are  three  exchange  current  densities,  one  from  the  reactions  at 
the  anode-electrolyte  interface,  another  from  the  reactions  at  the 
cathode-electrolyte  interface  and  the  third  from  the  reactions  from 
the  cathode-pore  interface. 

For  closed-circuit  operation,  we  continue  to  use  the  concept 
of  forward  and  reverse  reactions  and  their  corresponding  current 
densities.  The  net  current  density  across  the  cell  comes  from  the  dif¬ 
ference  between  the  forward  and  reverse  current  densities,  each  of 
which  is  described  as  the  product  of  an  attempt  current  density  and 
a  reaction  success  probability. 

In  our  model,  we  treat  all  the  reactions  at  interfaces  as  collision 
events  which  can  be  described  by  classical  (Boltzmann)  statistics. 
The  forward  current  density  at  the  anode-electrolyte  interface  can 
be  described  as 

iaf  =  ae01  /2  cosh  a,  and  a  =  {Ua  -  qVa)/2kT 
Here  a  is  forward  attempt  current  density  and  e“/2cosho'  is 
reaction  success  probability  for  the  reaction,  H2^2H+  +  2e_,  at 
anode-electrolyte  interface.  Ua  is  the  enthalpy  released  when 
the  reaction  occurs,  q  is  the  charge  transfer  per  reaction  which 
is  3.2  x  lO-19  C,  and  Va  is  the  potential  difference  across  the 
anode-electrolyte  interface.  Thus,  ( Ua-qVa )  is  the  net  energy 
required  for  the  reaction  to  occur.  The  expression  for  a  is 

1  /  JT  \  1/2 

a  =  2nH2  (  )  faQ(  1  -  %e)D  -  a3elenHMe/3)2,  (9) 

where  mH2  is  the  H2  molecule  mass,  and  ( kT/m H2f^2  is  the 
H2  molecule  average  velocity  component  (through  the  anode 
pore)  directed  toward  the  triple  phase  boundary  (TPB).  (1  - 
veie){  1  -  a3enH,de/3)2  is  the  probability  that  an  O2-  ion  without 
a  proton  attached  to  it  at  the  anode  TPB  so  that  the  H2  dissociation 
reaction  can  occur  in  which  vde  is  the  probability  of  oxygen  ion 
vacancy  in  electrolyte  material,  a3  e/ 3  is  the  unit  volume  of  oxygen 
ion  in  unit  lattice  of  electrolyte  material  and  nH)e/e  is  proton  con¬ 
centration  in  the  electrolyte,  and  fa  is  the  fraction  of  oxygen  sites 
from  electrolyte  that  sit  on  the  TPB  and  are  available  for  the  reac¬ 
tion  at  the  anode-electrolyte  interface.  The  detailed  derivation  and 
discussion  of  our  V(i)  model  will  be  published  elsewhere. 

By  setting  nn2  =  Hh2,  as  under  the  anode  saturation  condition,  Eq. 
(8c)  can  be  rewritten  as 
nH2,p,  ~  CH2,asT3  —  nN2,p,as(exp(CN2,as'r3)  —  1) 

-  nH2O,p,as(exp(CH2o,asT3)-l)-nH2,as  =  0  (10) 

All  constants  in  Eq.  (10)  are  positive  and  can  be  solved  numeri¬ 
cally  for  r. 

The  constants  in  Eq.  (10)  have  been  defined  except  for  the 
boundary  conditions  nH2,p,as,  nN2,P,as  and  nH2o,P,as  for  gas  concen¬ 
trations  in  the  plenum  at  the  outer  anode  surface  and  nH2,as  for  the 
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H2  concentration  at  the  anode-electrolyte  interface.  To  find  nH2,p,as 
we  consider  the  total  flow  in  and  out  of  the  anode  plenum.  For 
all  the  experiments,  the  total  gas  inflow  rate  before  going  through 
the  water  bubbler  is  200  ml  min-1.  The  total  metered  inflow  rate 
of  200  ml  min-1  at  temperature  21  °C  and  at  Bozeman,  Montana 
atmosphere  (4900  ft.  above  sea  level,  Peozeman  =  8.534  x  104  N  m-2) 
is,  from  the  ideal  gas  law, 

^Bozeman  S0  =Jl  =  7.012  x  1019  molec  S_1 

(11) 


The  metered  inflow  H2  is  found  by  multiplying  the  total  metered 
inflow  rate  by  the  partial  pressure  fraction,  and  N2  is  the  bal¬ 
ance  of  the  total  inflow.  Because  H20  is  obtained  by  flowing 
the  inflow  gas  through  a  water  bubbler,  the  flow  rate  of  H20  is 
the  inflow  rate  multiplied  by  saturated  water  partial  pressure  at 
room  temperature  (3%).  For  example,  a  20%  FI2  input  gas  flow 
has  metered  flow  rates:  ju2,m  =  Ph2Ji  =  1-404  x  1019  molec  s-1, 
jN2,m  =  Pn2Ji  =  5.618  X  1019  molec  s_1  andjH2o,m  =  Ph2o,scJi  = 
2.107  xlO18  molec  s-1. 

To  find  the  net  flow  of  each  gas  into  the  plenum,  we  must 
consider  the  fuel  gas,  FI2,  outflow  into  the  anode.  The  inflow  rate 
of  H2  into  anode  is  j’h2  =  - iS/q ,  where  S  is  the  active  area  and 
q  =  3.2  x  1 0-19  C.  When  the  cell  is  running  at  saturation  current,  the 
inflow  rate  into  the  anode  is  designated  as  Jh2, as- To  summarize  the 
analysis  for  net  inflows  into  the  plenum, 

iS 

jH2,net  =  P\\2J\  ~  Ju2,net  =PNjl,  ^2) 

Jh20,  net  =  Pu20,sad\  and  j  total, net  =  JH2,net  +  jN2,net  +Jh20,  net 

Our  test  system  is  open  to  space  and  we  assume  that  the  gases  in 
the  plenum  are  well  mixed  and  have  the  same  mole  fraction  every¬ 
where.  Accordingly,  in  the  plenum  the  gas  concentrations  under 
anode  saturation  conditions  are 


Pbo 


fttotal,p 
ftN  2,p,as  — 
t*H20  ,p,as  — 


kJ 

]N2,net 
J  total, net,  as 

j\i20,net 
J  total, net, as 


nH2,p, 


Ju2,net,as 


J  total,  net,  as 
fttotal,p  and 

ft-total,p 


fttotal,p 


(13) 


The  remaining  parameter  to  determine  is  nH2,as>  the  H2  con¬ 
centration  at  the  anode-electrolyte  interface  under  the  anode 
saturation  condition.  From  Eq.  (9),  nH2  remains  proportional  to  a  as  i 
changes,  because  we  assume  the  other  parameters  are  independent 
of  i.  For  open-circuit  (i  =  0)  conditions,  Hh2  at  the  interface  equals 
nH2,o  =  Ph2  ft  total,  p-  This  known  value  for  open-circuit  nn2  inserted 
into  Eq.  (9)  provides  a  known  value  a0  =  Ph2  aoi  for  the  open-circuit 
value  of  a,  where  a0i  is  the  open-circuit  value  for  pH2  =  1  and  is 
51 0.60  and  486.22  A  cm-2  for  700  and  800  °C,  respectively.  For  the 
anode  saturation  condition,  a  =  ias,  and  because  Hh2  is  proportional 
to  a,  we  have 


nH  2,  as  = 


Ph2  ft  total, p^as  _  Ktotal,plas 


a0 


P  H2Qoi 


CZqi 


(14) 


Inserting  the  nn2,as  value  and  the  nH2,p,as  value  into  Eq.  (10) 
yields  the  following  equation  for  r3  in  terms  of  known  parameters 


jH2,net,as 
J  total,  net,  as 


fttotal,p  ~  Q-I2,as^3  —  nN2,P,as(exp(cN2?asr3)  —  1) 


-nH20,p,as(exp(CH20,as'!3)  “  1) 


ft  total, p^as 


=  0 


4.  Results  and  discussion 

For  the  purpose  of  tortuosity  investigation,  a  series  of  differ¬ 
ent  thickness  anode-supported  Ba(Ceo.8Yo.2)03_«$  FI-SOFCs,  from 
1.85  to  1mm,  were  made  and  tested.  Fig.  2  shows  the  cross- 
section  microstructure  of  one  of  the  tested  Ba(Ceo.8Yo.2)03_«$ 
anode-supported  cells  from  SEM.  The  dense  electrolyte  is  ~20  [xm 
in  thickness  and  adheres  to  a  ~35  |xm  anode  interlayer  very  well. 
The  thickness  of  porous  cathode  is  ~45  |xm.  The  mean  pore  radius  of 
the  supporting  anodes  estimated  quantitatively  from  SEM  micro¬ 
graphs  was  ~1.41  |xm.  The  porosity  of  the  supporting  anode  was 
~34%,  which  was  measured  using  Archimedes’  method  described 
in  the  experimental  procedures.  This  porosity  only  refers  to  open 
pores  of  the  supporting  anode  since  closed  pores  do  not  provide 
channels  for  gas  diffusion. 

The  electrochemical  performances  of  a  Ba(Ce0.8Yo.2  P3-S  anode- 
supported  H-SOFC  with  a  1.5  mm-thick  anode  are  shown  in  Fig.  3. 
The  tests  were  done  at  temperature  700  °C,  Fig.  3(a),  and  800 °C, 
Fig.  3(b),  under  different  FI2  partial  pressures  while  air  was  used 
as  oxidant  at  the  cathode.  The  near  linear  l-V  curves  of  the  cell 
performance  at  high  H2  partial  pressures,  >40%  H2,  indicate  the 
domination  of  ohmic  polarization  of  the  cell  outputs.  For  low  FI2 
partial  pressures,  less  than  30%,  the  convex  l-V  curvatures  indicates 
the  activation  and  ohmic  polarization  loss  of  the  cell  are  small  com¬ 
pared  to  the  concentration  polarization,  especially  when  the  cell 
was  running  at  high  current  density.  The  concentration  polariza¬ 
tion  effect  from  changing  FI2  partial  pressures  also  shows  up  with 
the  decreasing  of  open-circuit  voltage  when  H2  partial  pressure 
decreases,  Fig.  4.  When  comparing  the  open-circuit  voltage  (OCV) 
of  H2/air  gas  inputs,  our  OCV  is  lower  than  the  theoretical  value, 
by  about  0.1  V  at  700  °C  and  0.122  V  at  800  °C.  This  is  because  of 
the  redox  reaction  of  Ce3+/Ce4+  in  Ba(Ce0.8Yo.2)03_5  under  reduc¬ 
ing  atmosphere,  which  increases  the  electronic  conductivity  of  the 
electrolyte  and  results  in  lower  OCV.  The  electronic  conductivity  of 
Ba(Ce0.8Yo.2  )03_5  increases  with  increasing  operating  temperature, 
which  explains  the  bigger  voltage  difference  of  OCV  between  the 
theoretical  EMF  and  measured  OCV  at  800  °C.  Similar  results  were 
reported  when  using  lanthanum-doped  ceria  and  zirconium-doped 
BCY  as  electrolyte  for  SOFCs  [14-19]. 

When  examining  the  l-V  curves  from  Ba(Ce0.8Yo.2)03_5  anode- 
supported  cells,  we  notice  the  current  density  increases  rapidly  at 
low  voltage,  e.g.  at  voltage  lower  than  0.05  V  for  1 0%  FI2  partial  pres¬ 
sure  in  Figs.  3(b)  and  Fig.  5(a).  The  same  feature  is  reproducible  and 
was  also  seen  on  the  anode  thickness  of  1.34  and  1.84  mm  but  not 


floi 


(15) 


Fig.  2.  Microstructure  of  Ba(Ceo.8Yo.2)03_5  anode  supported  H-SOFC  from  SEM 
image. 
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■  100%H2 


Current  density  (mA/cm2) 

■  ioo%h2 


Current  density  (mA/cm2) 

Fig.  3.  Hydrogen  partial  pressure  dependent  electrochemical  performance  of 
Ba(Ceo.8Yo.2)03_5  anode  supported  H-SOFC  with  1.5  mm  anode  thickness  at  (a) 
700 °C  and  (b)  800 °C. 


on  1  mm  cells.  It  seems  the  “tail”  gets  bigger  when  the  support¬ 
ing  anode  gets  thicker.  The  cell  with  1.5  mm  anode  thickness  was 
specifically  run  with  20%  H2  partial  pressure  under  0.1  V  for  more 
than  13  h  at  800  °C  before  the  test  system  was  shut  down.  Since 
the  cell  can  be  run  for  a  long  period,  it  excludes  the  possibility  that 
the  electrolyte  was  giving  up  its  oxygen  ion  for  the  SOFC  reaction, 
i.e.  the  electrolyte  was  not  decomposing  and  cannot  provide  the 
cell  an  extra  oxygen  source.  Conventionally,  we  expect  a  straight 
down  I-V  curve  when  concentration  polarization  dominates,  using 


Fig.  4.  Open-circuit  voltage  of  Ba(Ce0.8Yo.2)03_5  anode  supported  H-SOFC  with 
1.5  mm  anode  thickness. 


up  all  of  the  fuel  gas  at  the  anode-electrolyte  interface,  and  giving 
a  saturation  current  density.  This  behavior  indicates  a  mechanism 
faster  than  gas  diffusion  dominated  supply  of  fuel  gas  to  the  TPB. 
Our  speculation  on  this  mechanism  is  surface  diffusion  of  protons 
along  Ni  surface  to  the  TPB.  For  FI-SOFC,  protons  diffusing  along 
Ni  surface  allows  them  to  incorporate  into  electrolyte  at  the  TPB 
and  provides  enough  protons  for  the  reaction  on  the  cathode  side, 
i.e.  2H+  +  O2-  H20.  The  same  mechanism  was  not  seen  on  our  O- 

SOFC,  8  mol%  yttria-stabilized  Zr02  cells,  because  two  protons  must 
meet  an  O2-  at  the  TPB  at  the  same  time  for  the  reaction  to  occur  on 
the  anode  side.  However,  Williford  et  al.  and  Lee  et  al.  point  out  that 
surface  diffusion  of  protons  along  Ni  surface  occurs  only  over  small 
distance  and  lower  speed  when  compared  to  bulk  gas  diffusion 
[4,8].  Therefore,  the  reason  for  the  “tail”  behavior  remains  unclear. 

Due  to  the  fact  that  a  was  obtained  from  the  modeling  of  imping¬ 
ing  rate  of  H2  onto  the  triple  phase  boundary  (TPB),  the  amount  of 
Hh2  present  at  the  anode-electrolyte  interface  determines  the  sat¬ 
uration  current  density  when  V  goes  to  zero.  Therefore,  we  took 
the  current  density  as  ias  when  the  voltage  drops  to  zero  in  the 
I-V  curve  measurements.  Values  of  ias  for  various  hydrogen  par¬ 
tial  pressures  appear  in  Table  2.  Fig.  5  shows  the  results  from  our 
tortuosity  calculations  and  the  electrochemical  performances  of  a 
Ba(Ceo.8Yo.2)03-<$  anode-supported  H-SOFC  with  an  anode  thick¬ 
ness  1 .84  mm.  The  shown  data  is  for  tests  at  800  °C  under  different 
H2  partial  pressures.  The  tortuosity  results  range  from  2.94  to  1.50, 
but  are  limited  to  a  value  around  1.50  when  H2  concentrations  are 
below  20%.  As  shown  in  Fig.  5(a),  the  cell  performances  did  not 
encounter  serious  concentration  polarization  at  H2  partial  pressure 
above  50%.  For  H2  partial  pressure  above  50%,  the  main  loss  was 


Table  2 

List  of  saturation  current  densities  for  different  anode  thicknesses  which  were  tested  under  different  H2  partial  pressures  at  (a)  800  °C  and  (b)  700  °C. 


Ph2 

ias  (A cm-2)  Anode  1.84mm 

ias  (A  cnrr2)  Anode  1.50  mm 

ias  (A  cm-2 )  Anode  1 .34  mm 

ias  (A  cm-2)  Anode  1.00mm 

(a) 

10% 

0.506 

0.522 

0.709 

0.854 

15% 

0.782 

0.816 

0.955 

1.174 

20% 

1.055 

1.091 

1.193 

1.410 

30% 

1.489 

1.537 

1.493 

1.686 

(b) 

10% 

0.543 

0.517 

0.619 

0.757 

15% 

0.721 

0.782 

0.802 

1.058 

20% 

0.846 

1.015 

0.950 

1.294 

30% 

1.005 

1.338 

1.117 

1.590 

698 


C.-L.  Tsai,  V.H.  Schmidt  /  Journal  of  Power  Sources  196  (201 1 )  692-699 


•  100%  H2 


(b) 


no 

100 

90 


cell  tested  temperature  800  °C 
anode  porosity  34% 
average  pore  radius  1.41  pm 
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Fig.  5.  Ba(Ce0.8Yo.2)03_5  anode-supported  H-SOFC  (a)  the  electrochemical  perfor¬ 
mance  at  800  °C  under  different  H2  partial  pressures  and  (b)  calculated  tortuosities 
for  its  anode. 


from  ohmic  polarization  so  that  the  calculated  tortuosity  using  ias 
at  voltage  drop  to  zero  should  not  be  considered.  Note  that  when 
concentration  polarization  is  severe,  the  tortuosity  values  approach 
the  same  value.  Therefore,  the  tortuosity  calculated  from  lowest  H2 
concentrations  should  be  a  more  accurate  number.  Fig.  5(b)  also 
shows  the  calculated  tortuosity  for  H2  partial  pressure  under  40% 
when  the  “tail”  effect  is  ignored,  using  the  intersect  value  between 
solid  lines  and  current  density  axis  in  Fig.  5(a).  The  tortuosity  is 
then  increasing  from  ~1.50  to  ~1.60  at  low  FI2  concentration. 

Ba(Ceo.8Yo.2p3-5  anode-supported  cells  with  different  anode 
thickness,  but  the  same  porosity  and  pore  radius,  were  tested  and 
the  tortuosities  of  their  anodes  were  also  calculated.  The  results 
are  shown  in  Fig.  6.  The  tortuosities  obtained  from  10%  H2  con¬ 
centration  all  fall  into  the  same  range,  1.55  ±  0.1,  and  are  thickness 
independent,  as  they  should  be,  for  the  tests  both  at  800  and  700  °C. 
Compared  to  700  °C,  the  tortuosity  values  approach  their  final  value 
at  higher  FI2  concentration  when  cells  were  tested  at  800  °C.  This  is 
because  of  the  higher  fuel  utilization  of  cells  at  800  °C  that  makes 
the  cells  encounter  concentration  polarization  at  higher  FI2  partial 
pressure. 

We  compare  the  result  presented  here,  1.55  ±  0.1  with  our  pre¬ 
vious  results  from  analyzing  Jiang-Virkar’s  experiment  which  was 
2.3  ±  0.6  [6,14].  If  we  take  into  account  that  the  enhancement  factor 
should  have  been  r2/0,  then  the  new  number  in  Ref.  [6]  is  1 .74  ±  0.3. 
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Fig.  6.  Calculated  tortuosities  for  different  thickness  anode  cells  tested  at  (a)  800  °C 
and  (b)  700  °C. 


The  tortuosity  value  presented  here  is  a  little  smaller  than  the  previ¬ 
ous  average  value  but  within  the  error.  This  results  from  our  bigger 
pore  radius,  1 .41  p,m,  when  compared  to  Jiang-Virkar’s  0.5  p,m  pore 
radius,  even  though  our  porosity  is  lower,  34%,  compared  to  their 
54%  porosity.  When  looking  at  Jiang-Virkar’s  20%FI2-80%N2-FI20 
ternary  system,  its  saturation  current  ias  is  0.89  A  cm-2.  The  test  in 
Fig.  4  with  20%  H2  partial  pressure  delivered  a  saturation  current 
1 .055  A  cm-2  or  0.877  A  cm-2  when  the  “tail”  is  not  considered.  For 
our  much  lower  porosity  but  delivering  similar  amount  of  FI2  fuel 
from  plenum  to  anode-electrolyte  interface,  it  makes  sense  that  our 
tortuosity  is  lower,  i.e.  the  pores  are  straighter  from  plenum-anode 
interface  to  anode-electrolyte  interface.  Furthermore,  the  bigger 
pore  size  reduces  the  collision  rate  of  gas  molecules  with  the 
wall,  thereby  reducing  the  effect  from  Knudsen  diffusion,  result¬ 
ing  in  lower  tortuosity  of  our  supporting  anode.  This  value  is 
also  similar  to  the  results  which  were  obtained  from  FIB-SEM,  X- 
ray  computed  tomography  and  gas  counter  diffusion,  r  =  1.33-4.0 
[3-5]. 

Plots  of  gas  concentrations  across  the  anode  for  six  different  FI2 
partial  pressures  under  anode  limiting  currents  at  700  °C  are  shown 
in  Fig.  7.  The  total  concentrations  only  decrease  slightly  because 
of  the  compensation  from  N2  concentration  increasing.  The  expo¬ 
nential  effects  on  the  gas  concentrations  appear  insignificant.  They 
appear  more  like  linear  changes,  due  to  the  relatively  “thin”  sup¬ 
porting  anode.  To  see  the  exponential  effects,  it  is  necessary  to  have 
a  supporting  anode  thicker  than  1  cm  and  high  FI2  partial  pressure. 
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Fig.  7.  Plots  of  gas  concentration  vs.  position  in  anode  for  1 .84  mm  anode  thickness  under  anode-limiting  current  conditions  at  700  °C. 


5.  Conclusions 

A  series  of  different  anode  thickness  Ba(Ce0.8Y0.2)O3_5  anode- 
supported  H-SOFC  were  made  and  tested.  By  changing  H2  partial 
pressures,  the  anodic  concentration  polarization  effect  can  be 
observed.  For  low  FI2  partial  pressures,  saturation  current  densi¬ 
ties  under  concentration  polarization  were  obtained  and  used  for 
tortuosity  calculation.  The  result  of  tortuosity  calculation  for  the 
supporting  anode  is  1 .55  ±  0.1 ,  which  is  temperature  and  thickness 
independent,  as  it  should  be  since  tortuosity  is  a  geometric  factor 
of  the  anode  structure. 

The  accurate  tests  bring  down  the  error  between  cells  to  ±  0.1 
of  the  total  value.  The  value  of  1.55  ±0.1  is  a  little  smaller  than 
the  average  number,  which  we  found  previously  from  analyzing 
the  Jiang-Virkar  results,  1.74  ±0.3.  Flowever,  the  lower  poros¬ 
ity  and  bigger  pore  size  of  our  supporting  anode  are  consistent 
with  a  lower  tortuosity  to  deliver  the  measured  similar  amount 
of  H2  gas  from  the  plenum  to  the  anode-electrolyte  interface.  Our 
result  is  also  consistent  with  the  values  found  by  different  tech¬ 
nologies  for  modern  anode  material  with  typical  porosity  >30%. 
Therefore,  we  can  say  that  our  value  is  in  the  physically  reasonable 
range.  The  developed  equations  also  give  information  for  pressure 
drops  across  the  electrodes  and  the  gas  concentrations  at  the  H- 
SOFC  electrode-electrolyte  interface  which  provides  the  required 
tools  for  modeling  the  voltage-current  density  expression.  Such 
voltage-current  modeling  gives  information  for  new  SOFC  designs 
and  for  analyzing  performance  of  existing  designs  in  both  fuel  cell 
and  electrolysis  modes. 
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